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A  Generic  Modeling  Approach  to  Biomass 
Dynamics  of  Sagittaria  lati folia 
and  Spartina  a  Item  i flora 

by  Elly  P.  H.  Best,  William  A.  Boyd,  and  Kevin  P.  Kenow 


PURPOSE:  This  technical  note  describes  an  ecological  modeling  approach  that  can  be  used  to 
explore  relationships  between  species  of  emergent  aquatic  vegetation  communities  and  their 
environmental  conditions.  The  modeling  approach  was  used  to  evaluate  the  potential  persistence 
of  two  desired  —  and  quantitatively  important  —  rhizomatous  plant  species  under  various  cli¬ 
matological  conditions:  Sagittaria  latifolia,  common  in  freshwater  systems,  can  produce  tubers 
as  well  as  rhizomes;  Spartina  alterniflora  is  typical  for  coastal  marshes.  Both  species  are 
endemic  to  the  United  States. 

BACKGROUND:  Emergent  aquatic  vegetation  may  play  important  roles  in  aquatic  ecosystems. 
Functions  attributed  to  “desirable”  species  are:  stabilization  of  sediment  and  shores,  amelioration 
of  water  transparency,  regulation  of  nutrient  availability  in  the  water  column  and  service  as  a 
habitat  and  food  source  for  invertebrates,  fish  and  waterfowl.  Conversely,  effects  attributed  to 
“nuisance”  or  “invasive”  species  are:  excessive  biomass  production  that  interferes  with  human 
utilization  of  water  resources,  or  displacement  of  desirable  indigenous  communities.  Distribution 
and  abundance  of  emergent  vegetation  in  large  water  bodies  in  the  Mississippi  River  System 
(UMRS)  and  Coastal  Louisiana  (CL)  have  changed  for  decades  and  declined  in  recent  years.  In 
general,  variation  in  environmental  factors  such  as  water  depth,  temperature,  clarity,  current, 
wave  action,  and  substrate  characteristics  would  be  expected  to  affect  the  distribution  and 
production  of  emergent  macrophytes  (Gosselink  and  Turner  1978).  Sagittaria  latifolia 
(Broadleaf  arrowhead)  is  a  desired  and  dominant  species  in  the  UMRS,  where  it  provides  a 
significant  annual  autochthonous  input  (Eckblad  et  al.  1977).  In  this  river  system,  changes 
resulting  from  the  man-made  modification  of  the  hydrologic  cycle  include  installing  a  system  of 
dams  in  the  1930s,  navigation  pools  with  artificially-maintained  high  water  levels,  island  loss 
due  to  erosion,  and  increased  sedimentation  (Bellrose  et  al.  1979;  Eckblad  et  al.  1977). 
Environmental  changes  (such  as  increased  water  level  and  turbidity)  resulting  from  the  operation 
of  this  navigation  system  for  barge  transportation  of  bulk  commodities,  have  been  listed  as 
contributing  to  the  decline  of  S.  latifolia.  The  latter  statement  is  confirmed  by  the  fact  that 
experimental  decreases  in  water  level  during  the  summer  growth  season  in  UMRS  Pool  5  led  to 
the  increased  abundance  and  distribution  of  emergent  vegetation  (Kenow  et  al.  2007).  Spartina 
alterniflora  (Smooth  cordgrass)  is  a  desired  species  in  CL,  where  it  dominates  large  portions  of 
the  salt  marshes  because  of  its  high  primary  production  (Kirby  and  Gosselink  1976;  Darby  and 
Turner  2008).  Its  spatial  distribution  is  limited  to  the  coastal  areas  along  the  entire  Atlantic  and 
the  Southern  Pacific  seaboards  of  the  United  States.  Large-scale  diebacks  have  occurred  and 
have  been  attributed  to  multiple  causes:  pennanent  dieback  to  prolonged  flooding  of  the 
subsiding  marsh  surface  (Webb  et  al.  1995);  temporary  dieback  to  drought  from  which  the 
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vegetation  could  recover  rapidly  by  regrowth  from  rhizomes  or  more  slowly  in  the  absence  of  the 
rhizomes  by  seedling  recruitment  in  the  opened  areas  (Edwards  et  al.  2005). 

Besides  light  (as  affected  by  water  level  and  turbidity),  nitrogen  (N)  and  phosphorus  (P)  are  gen¬ 
erally  believed  to  be  the  most  important  limiting  factors  in  aquatic  systems  (Hutchinson,  1975). 
Relationships  between  biomass  nutrient  concentrations  and  nutrient  limitation  are  complex. 
Biomass  nutrient  concentrations  tend  to  be  positively  correlated  with  nutrient  supply  when  all 
other  resources  are  sufficiently  available  (Guesewell  and  Koersehnan  2002).  A  low  concentra¬ 
tion  of  N  in  plant  biomass  should  reflect  a  low  availability  of  N  to  this  plant  and,  therefore,  indi¬ 
cate  that  additional  supply  of  N  would  increase  the  plants’  biomass  production.  By  definition, 
this  means  that  N  is  limiting  (Vitousek  and  Howarth  1991).  If  two  or  more  nutrients  (e.g.,  N  and 
P)  are  in  short  supply,  their  availability  relative  to  each  other  is  likely  to  determine  which  of 
them  is  limiting.  Thus,  the  molar  ratio  of  N:P,  rather  than  the  individual  concentrations,  should 
indicate  limitation  (Koersehnan  and  Meuleman  1996),  with  tissue  N:P  ratios  less  than  14  indic¬ 
ative  of  N-limited  growth  in  terrestrial  plants  (Aerts  and  Chapin  2000).  Nutrient-related  growth 
limitation  of  emergent  plants  in  natural  systems  has  been  reported,  but  N:P  ratios  have  not 
usually  been  determined.  In  S.  latifolia,  N  may  limit  growth  under  natural  conditions  since 
results  of  a  short-term  pot  experiment  indicated  that  this  plant  depleted  the  exchangeable  N  in 
natural  sediments  within  four  weeks,  but  left  substantial  P  levels  (Barko  et  al.  1988;  Barko  and 
Smart  1983).  In  S.  alterniflora,  it  was  suggested  that  N  also  limits  growth  under  natural  condi¬ 
tions,  based  on  the  vegetation  response  to  N  fertilization  in  the  field  (Valiela  and  Teal  1974; 
Gallagher  1975).  In  contrast,  the  average  N:P  ratio  of  16:1  in  aboveground  biomass  and  of  37:1 
in  belowground  biomass  suggested  P  limitation  in  belowground  biomass  at  this  CL  site  (Darby 
and  Turner  2008).  However,  growth  limitation  by  nutrient  availability  may  be  even  more  com¬ 
plicated  in  plants  exposed  to  different  salinity  levels,  as  illustrated  by  the  results  of  a  pot  experi¬ 
ment  in  which  the  critical  N:P  ratios  were  determined  in  S.  alterniflora.  In  these  plants,  growth 
was  limited  by  N  with  tissue  N:P  ratios  <13  and  aboveground  biomass  was  correlated  with 
interstitial  sediment-N  concentration,  but  growth  rate  was  affected  by  salinity  (Smart  and  Barko 
1980).  The  relationships  between  the  potential  persistence  of  both  plant  species,  nutrient  limita¬ 
tion,  and  salinity  have  not  yet  been  unequivocally  elucidated;  consequently,  placeholders  have 
been  included  in  the  model  in  anticipation  of  supporting  data. 

Simulation  models  that  include  descriptions  of  aquatic  vegetation  responses  to  changes  in  physi¬ 
cal  and  chemical  conditions  in  various  climates  can  be  valuable  tools  for  water  resource  manag¬ 
ers.  These  models  can  be  used  to  evaluate  key  environmental  conditions  in  which  the  vegetation 
would  persist  or  produce  excessive  biomass,  with  ensuing  consequences  for  the  systems  in  which 
they  grow.  Additionally,  the  models  may  provide  insight  as  to  how  the  vegetation  would  be 
affected  by  different  management  scenarios  (Carr  et  al.  1997;  Best  et  al.  2001;  Karunaratne  and 
Asaeda  2002;  Asaeda  et  al.  2008).  In  this  paper,  a  dynamic  simulation  modelling  approach  to 
emergent  plant  biomass  formation  is  summarized,  with  light  and  temperature  as  driving  vari¬ 
ables,  and  including  descriptions  of  plant  responses  to  human  influences  such  as  management 
measures  resulting  in  changes  in  turbidity,  mechanical  harvesting,  grazing,  and  flooding.  Cali¬ 
bration  of  plant  responses  to  current  velocity  and  nutrient  limitation  will  be  added  later  on  when 
calibration  values  become  available.  This  modelling  approach  was  applied  to  submersed  aquatic 
vegetation  (SAV)  also  (Best  and  Boyd  2008).  The  approach  is  mathematically  similar  to  those 
followed  in  other  models  for  emergent  vegetation,  such  as  those  developed  for  Phragmites  aus¬ 
tralis  (Ondok  1973;  Asaeda  and  Karunaratne  2000;  Asaeda  et  al.  2008)  and  for  S.  alterniflora 


2 


ERDC  TN-SWWRP-1 1-1 
January  201 1 


(Morris  1989;  Dai  and  Wiegert  1996a).  The  approach  describes  plant  morphology  and  biomass 
formation  in  relative  detail,  but  it  differs  from  them  in  that  it  relates  ecophysiological  processes 
to  developmental  cycle,  using  the  model  to  simulate  plant  communities  in  different  climates.  In 
this  construct,  aerial  shoots  absorb  CO2  from  air,  and  submerged  shoots  absorb  it  from  water 
where  CO2  availability  is  assumed  to  be  typical  for  hard  water  with  an  alkalinity  between  0  and 
300  mg  L'1  and  a  circumneutral  pH;  effects  of  changes  in  CO2  availability  are  not  included.  The 
model  species  are  S.  latifolia  and  S.  alterniflora',  both  plants  have  similarities  in  growth  strategy, 
but  are  significantly  different  in  morphology  and  physiology.  The  model  has  been  calibrated, 
tested  for  sensitivity  and  validated  against  field  data  for  both  species.  Both  species  have  the 
capacity  to  persist  in  eutrophic,  shallow  water  bodies  with  fluctuating  water  levels;  both  are 
capable  of  forming  substantial  above-and  belowground  biomass,  thereby  functioning  as  impor¬ 
tant  marsh-characteristic  elements.  Important  physiological  differences  are  that  S.  latifolia  fixes 
carbon  via  the  C3  photosynthetic  pathway  in  contrast  to  S.  Alterniflora,  which  uses  the  C4  path¬ 
way;  S.  latifolia  has  a  higher  potential  photosynthetic  rate  at  light  saturation  for  aerial  shoots,  and 
a  higher  species-characteristic  light  extinction  coefficient  than  S.  alterniflora',  S.  latifolia  can 
form  tubers  (i.e.,  organs),  enabling  survival  during  adverse  conditions  such  as  drought  and  cold; 
S.  latifolia  is  sensitive  to  increased  salinity.  These  are  characteristics  which  make  S.  latifolia  a 
species  that  may  predominate  in  freshwater  and  coastal  marshes  with  a  relatively  low  salinity, 
whereas  S.  alterniflora  grows  well  in  coastal  marshes  of  higher  salinity.  With  both  species  hav¬ 
ing  different  strategies  to  survive  adverse  conditions,  changes  in  spatial  distribution  and 
replacement  of  one  species  by  the  other  (the  latter  in  coastal  areas  only)  over  a  period  of  one  to 
several  years  can  be  expected.  Consequently,  the  described  ecological  model  provides  an  ideal 
means  to  investigate  the  effects  of  relatively  short-term  changes  in  environmental  conditions  on 
the  potential  persistence  of  these  two  emergent  species  in  shallow  water  bodies  as  part  of  resto¬ 
ration  plans,  provided  detailed  information  on  environmental  conditions  is  available. 

Besides  model  calibration  and  validation,  two  other  aspects  of  the  relationship  between  impor¬ 
tant  representatives  of  rhizomatous  plant  species  and  environmental  conditions  were  investigated 
in  the  present  study.  A  dynamic  ecological  modeling  approach  was  used:  (i)  persistence  at  vari¬ 
ous  flood  and  drought  conditions,  and  (ii)  persistence  under  more  southern  climatological  condi¬ 
tions  than  at  the  calibration  site. 

ECOLOGICAL  MODELING  APPROACH:  This  ecological  model  type  simulates  the  carbon 
flow  mass  balance  of  typical  emergent  vegetation  on  a  1-m2  sediment/soil  with  an  overlying 
water  column  (Figure  1).  Growth  is  considered  as  the  plant  dry  matter  accumulation,  including 
rhizomes,  and,  if  present,  subterranean  tubers,  in  an  environment  where  N  and  P  may  be  limiting 
under  the  prevailing  weather  conditions.  At  least  one  plant  cohort  waxes  and  wanes  per  season  in 
different  climatological  conditions,  varying  from  temperate  to  tropical.  The  rate  of  dry  matter 
accumulation  is  a  function  of  irradiance,  temperature,  CO2  availability  and  plant  characteristics. 
The  rate  of  CO2  assimilation  (photosynthesis)  of  the  plant  community  depends  on  the  radiant 
energy  absorbed  by  the  canopy.  The  daily  rate  of  gross  CO2  assimilation  of  the  community  is 
calculated  from  the  absorbed  radiation,  the  photosynthetic  characteristics  of  leaves  and  the  CO2 
availability.  Calculations  are  executed  in  a  set  of  subroutines  added  to  the  model.  Part  of  the  car¬ 
bohydrates  produced  is  used  to  maintain  the  existing  biomass.  The  remaining  carbohydrates  are 
converted  into  structural  dry  matter  (plant  organs).  In  the  conversion  process,  part  of  the  weight 
is  lost  in  respiration.  The  dry  matter  produced  is  partitioned  among  the  various  plant  organs 
using  partitioning  factors,  defined  as  a  function  of  the  phenological  cycle  of  the  community.  The 
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dry  weights  (DW)  of  the  plant  organs  are  obtained  by  integration  of  their  growth  rates  over  time. 
The  plant  may  over-winter  through  rhizomes  and/or  tubers  in  the  sediment  without  or  with  plant 
biomass  present.  Rhizomes  may  persist  when  a  critical  rhizome  mass  is  maintained.  Tubers  are 
depleted  and  disintegrate  in  the  summer  following  the  season  in  which  they  were  formed.  All 
calculations  are  perfonned  on  an  m  basis.  Since  environmental  factors  and  plant  growth  char¬ 
acteristics  vary  with  plant  height  and  water  depth,  in  the  model  the  growth-related  processes  of 
the  aboveground  plant  biomass  and  the  water  column  have  been  partitioned  in  0.1 0-m  depth  lay¬ 
ers.  A  relational  diagram  is  presented  in  Figure  2.  Seed  formation  has  not  been  included  in  the 
model,  because  its  role  in  maintaining  established  emergent  plant  communities  in  a  temperate 
climate  is  minimal.  Dispersal  and  colonization  of  new  habitats  are  recognized,  important  char¬ 
acteristics  of  emergent  plants.  The  latter  processes,  however,  are  better  described  using  other 
modelling  approaches  (based  on  logistic  regression  or  on  descriptions  of  population  dynamics 
varying  in  time  and  space),  as  described  by  Scheffer  (1991). 


Plant  processes 
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to  plant  organs 
•Photosynthetic 
activity  air  (C3-C4) 
•Photosynthetic 
activity  under  water 
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•N,  P  biogeochemistry 
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Figure  1 .  Schematic  generic  emergent  plant  growth  model. 


General  features  of  the  model  include  that  it: 

•  Is  operational  in  a  one-dimensional  (quasi  two-dimensional)  configuration 

•  Follows  a  state  variable  approach 

•  Provides  that  the  state  variable  selected  may  be  individually  activated  or  deactivated 

•  Performs  integration  using  the  Runge-Kutta  method 

•  Computes  photosynthesis  per  second  and  other  masses  per  day 
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TEMPERATURE  LIGHT  C02  AVAILABILITY  PHOTOSYNTHETIC 

EFFICIENCY 


Figure  2.  Relational  diagram  illustrating  the  following  model  processes  in  ARROW:  (1)  phonological 
cycle  and  development;  (2)  photosynthesis,  respiration,  and  biomass  formation;  and 
(3)  flowering,  translocation,  senescence  and  wintering  organs  (the  latter  process  in  grey 
background).  Rectangles  represent  quantities  (state  variables);  valve  symbols,  flows  (rate 
variables);  circles,  auxiliary  variables;  underlined  variables,  driving  and  other  external 
variables;  dashed  lines,  information  flow  (symbols  according  to  Forrester  1961). 
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•  Operates  as  a  stand-alone  version  fitted  in  a  FORTRAN  Simulation  Environment  (FSE) 
shell  (Van  Kraalingen  1995).  Provides  binary  and  ASCII  output  files,  and  graphics  that 
can  be  viewed  within  a  user-friendly  shell.  Coded  in  ANSI  Standard  FORTRAN  F77. 

Central  Features.  Central  features  of  the  model  are  the  (1)  link  between  the  species-charac¬ 
teristic  phenological  cycle,  physiological  processes  and  environmental  conditions  and  (2)  state 
variable  equation  detennining  instantaneous  gross  photosynthesis. 

Species-characteristic  Phenological  Cycle.  The  phenology  of  the  plant  community,  for 
which  the  development  phase  can  be  used  as  a  measure,  is  modelled  as  a  sequence  of  processes 
that  take  place  over  a  period  of  time,  punctuated  by  more  or  less  discrete  events.  The  develop¬ 
ment  phase  (DVS)  is  a  state  variable  in  the  models.  The  DVS  is  dimensionless  and  its  value 
increases  gradually  within  a  growing  season.  The  development  rate  (DVR)  has  the  dimension  d'1. 
The  multiple  of  rate  and  time  period  yields  an  increment  in  phase.  The  response  of  DVR  to  tem¬ 
perature  in  the  model  is  in  accordance  with  the  degree-day  hypothesis  (Thornley  and  Johnson 
1990).  Calibration,  according  to  this  hypothesis,  allows  use  of  the  model  for  the  same  plant  spe¬ 
cies  at  various  sites  with  different  climates  (temperature  regime).  The  relationships  between  the 
development  phase,  day-of-year,  and  3°C  day-degree  sum  for  a  temperate  climate  are  presented 
in  Table  1. 

Each  simulation  starts  at  the  first  Julian  day  (i.e.,  1  January,  when  the  DVS  has  the  value  of  0.0). 
For  S.  latifolia,  a  species  that  may  overwinter  with  rhizomes  and  tubers,  the  simulation  starts 
using  a  selected  rhizome  weight  and/or  tuber  bank  density/individual  tuber  weight  combination 
as  initial  values.  Initiation  of  growth  activity  occurs  by  sprouting  of  the  tubers,  or  sprouting  of  a 
fixed  number  of  plants  at  a  DVS  >  0.292.  Sprouts  of  the  first  plant  cohort  develop  through  remo¬ 
bilization  of  carbohydrates  until  the  tubers  or  rhizomes  are  depleted.  If  the  first  plant  cohort  does 
not  succeed  in  becoming  self-supporting  and  DVS  is  less  than  1.001,  a  second  cohort  sprouts 
from  the  tuber  bank  or  rhizomes.  For  S.  alterniflora,  a  species  that  overwinters  with  rhizomes, 
the  simulation  starts  using  a  selected  rhizome  weight.  The  DVS  values  of  the  phenological  proc¬ 
esses  in  S.  alterniflora  differ  from  those  in  S.  latifolia  (Table  1). 

Instantaneous  Gross  Photosynthesis  and  Biomass  Formation.  Light  availability  is  an 
important  factor  influencing  the  distribution  and  abundance  of  aquatic  plants.  For  emergent  leaves, 
light  attenuation  only  within  the  plant  canopy  occurs.  Leaves  submersed  in  water  may  have  a  small 
part  of  the  irradiance  reflected  by  the  water  surface,  and  further  attenuation  may  occur  by  water 
and  its  suspended  solids  and  by  the  plant  itself,  either  covered  by  epiphytes  or  not.  Emergent  leaves 
fix  carbon  with  a  higher  potential  photosynthetic  rate  at  light  saturation  than  submersed  leaves. 
Measured  daily,  total  irradiance  (wavelength  300-3000  nm)  is  used  as  input  in  the  model.  Only  half 
of  the  irradiance  reaching  the  water  surface  is  considered  to  be  photosynthetically  active  and  is, 
therefore,  used  as  a  base  for  the  calculation  of  CO2  assimilation.  Part  of  the  irradiance  (6  percent) 
can  be  reflected  by  the  water  surface.  The  subsurface  irradiance  can  be  attenuated  by  dissolved 
substances  and  particles  (in  mg  L'1)  within  the  water  column  resulting  in  a  site-  and  season-specific 
water  extinction  coefficient  (Equation  1).  The  remaining  radiation  may  be  further  reduced  by  epi¬ 
phyte  shading  (Equation  1.1).  The  vertical  profiles  of  the  radiation  within  the  plant  layers  are  cha¬ 
racterized  also.  The  absorbed  irradiance  for  each  horizontal  plant  layer  is  derived  from  these  pro¬ 
files  (Equation  1.2).  The  plant  light  extinction  coefficient,  K,  is  plant  species-characteristic  and 
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Table  1.  Relationship  between  plant  development  phase  (DVS),  day  of  year,  and 
3°C  day-degree  sum  in  a  temperate  climate3  (DVRVT=  0.015;  DVRRT=  0.040  for 
tuberforming  species;  DVRVT=  0.022;  DVRRT=  0.015  for  non-tuberforming 
species;  at  a  reference  temperature  of  30°C). 

Plant  developmental  phase 

S.  latifolia 

S.  alterniflora 

Description 

DVS  value 

Day 

number 

3°C  Day- 
degree  sum 

DVS  value 

Day 

number 

3°C  Day- 
degree  sum 

First  Julian  day  number  -> 
sprouting,  initiation  elongation, 
leaf  expansion  COHORT1 

0  0.291 

0  129 

1  341 

0  ->  0.375 

0  ^  42 

1  392 

Sprouting,  initiation  elongation, 
leaf  expansion  floral  initia¬ 
tion,  anthesis,  induction  of 
tuber  formation*  and  senes¬ 
cence  COHORT1 

0.292  ^  1.000 

130  -»  199 

342  ^  1682 

0.376  ^  1.000 

43  ^  99 

393  1080 

Floral  initiation,  anthesis, 
induction  of  tuber  formation* 
and  senescence  ->  transloca¬ 
tion,  tuber  formation*  and 
senescence  COHORT1 

1.001  1.630 

199  212 

1683  1955 

1.001  1.630 

100  154 

1081  -»  2171 

Translocation,  tuber  formation* 
and  senescence  ->  senesced 
COHORT1 

1.631  ^  2.000 

213  -»  321 

1956  2293 

1.631  ^  2.000 

155  183 

2172  ^  2827 

Sprouting,  initiation  elongation, 
leaf  expansion  floral  initia¬ 
tion,  anthesis,  induction  of 
tuber  formation*  and  senes¬ 
cence  COHORT2 

1.001  ^  1.630 

199  -»  212 

1683  1955 

1.001  -»  1.630 

100  154 

1081  -»  2171 

Floral  initiation,  anthesis, 
induction  of  tuber  formation* 
and  senescence  ->  transloca¬ 
tion,  tuber  formation*  and 
senescence  COHORT2 

1.631  2.000 

213  -»  321 

1956  2273 

1.631  2.000 

155  183 

2171  -»  2827 

Translocation,  tuber  formation* 
and  senescence  ->  senesced 
COHORT2 

2.001  2.570 

321  365 

2274  3409 

2.001  2.570 

184  ^  365 

2828  -»  6257 

Senesced  COFIORT  1  and  2 

2.570 

365 

3409 

2.570 

365 

6257 

a  Calibration  was:  for  S.  latifolia  on  field  data  on  biomass,  water  transparency  and  depth  from  Upper  Mississippi  River  Pool  9, 

Iowa,  1982  (Clark  and  Ciay  1985),  irradiance  from  La  Crosse,  Wisconsin,  1982,  water  temperature  from  Lansing,  Iowa;  for 

S.  alterniflora  on  field  data  on  biomass  from  Sapelo  Island,  Georgia  (Dai  and  Wiegert  1996a,b),  irradiance  and  air  temperature 
from  Brunswick,  Georgia,  1991,  and  fixed  water  transparency  and  depth.  *  Tuber  formation  only  in  S.  latifolia. 

assumed  to  be  constant  throughout  the  year.  The  incoming  irradiance  is  attenuated  by  the  shoots, 
part  of  which  is  absorbed  by  the  photosynthetic  plant  organs:  i.e.,  the  leaves.  Instantaneous  rates  of 
gross  assimilation  are  calculated  from  the  absorbed  light  energy  and  the  photosynthesis  light 
response  of  individual  shoots,  here  used  synonymously  to  leaves.  The  photosynthesis-light 
response  of  leaves  is  described  by  Equation  1.3.  In  the  photosynthesis-light  response  equation,  the 
value  of  potential  photosynthetic  activity  at  light  saturation  ( AMX)  is  characteristic  for  submersed 
plants,  and  the  AMX2  and  the  initial  light-use  efficiency  (EE)  typical  for  C3  plants  in  S.  latifolia 
and  for  C4  plants  in  S.  alterniflora.  AMX  and  AMX2  are  affected  by  temperature  via  a  fitted,  rela¬ 
tive  function,  AMTMPT,  accounting  for  the  measured  effect  of  daytime  temperature,  and  enabling 
the  calculation  of  the  actual  photosynthesis  rate  (AMAX).  AMAX  can  be  affected  by  tissue  N:P  ratio 
via  a  species-characteristic,  fitted,  relative,  function  NPREDF,  accounting  for  the  still-to-be- 
measured  effect  of  tissue  N:P  ratio  on  plant  biomass  production,  here  used  synonymously  for 
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photosynthesis.  AMAX  may  also  be  affected  by  current  velocity  via  a  species-characteristic,  fitted, 
relative,  function,  REDAM  1,  accounting  for  the  measured  effect  of  current  velocity  on  AMX. 
Senescence  may  affect  AMX.  Substituting  the  appropriate  value  for  the  absorbed  photosyntheti- 
cally  active  radiation  yields  the  assimilation  rate  for  each  specific  shoot  layer.  The  instantaneous 
rate  of  gross  assimilation  over  the  height  of  the  vegetation  is  calculated  by  relating  the  assimilation 
rate  per  layer  to  the  species-characteristic  biomass  distribution  and  by  subsequent  integration  of  all 
vegetation  layers.  The  daily  gross  assimilation  rate  is  calculated  by  using  the  Gaussian  integration 
method.  A  portion  of  the  carbohydrates  fonned  is  respired  in  maintenance  of  existing  plant  com¬ 
ponents  and  during  the  formation  of  new  plant  components  (i.e.  growth).  After  flowering,  down¬ 
ward  translocation  of  assimilates  start  filling  the  rhizomes;  tubers  may  be  induced  and  formed 
under  a  specific  combination  of  temperature  and  day  length,  and  senescence  sets  in.  The  S.  latifolia 
model  application  (ARROW)  was  calibrated  on  data  pertaining  to  a  S.  latifolia  vegetation  in  Upper 
Mississippi  River  System  Pool  9,  IA,  USA,  1982  (Clark  and  Clay  1985).  The  S.  alterniflora  model 
application  (CORDG)  was  calibrated  on  data  pertaining  to  a  S.  alterniflora  vegetation  on  Sapelo 
Island,  Georgia,  USA,  1991  (Dai  and  Wiegert  1996a,  b).  The  models  simulated  the  dynamics  of 
plant,  rhizome  and  tuber  biomass  and  tuber  numbers,  the  latter  for  S.  latifolia  only,  for  the  calibra¬ 
tion  and  validation  sites  well  over  a  period  of  one  to  five  years.  The  models  have  been  used  to 
simulate  plant,  rhizome  and  tuber  biomass  and  tuber  numbers  for  other  sites  with  temperate  and 
tropical  climates  as  well.  Key  model  equations  dealing  with  photosynthesis  and  nutrient  limitation 
are  provided  in  Appendix  A,  and  parameters,  variables,  and  constants  are  provided  in  Table  2. 
More  detailed  descriptions  of  the  equations  involved  and  model  applications  can  be  found  in  Best 
and  Boyd  (2007,  2008).  Executable  versions  of  the  models  are  available  at  http: //el. erdc. usace. 
army.mil/products.  cfin  ?  Topic  =model&  Tvpe=aquatic. 


IRZm  =  IRZ,.  x  e 


(-TLxL-KxSCj) 


(1) 


IABS,  = 


( IRZi  -  IRZ, , , )  x  SCI  x  K 
(. KxSCi+TLxL ) 


x(]  .0-  EPISIID) 


(1.1) 


1ABSL-  -  IABS,  x  FL 


FGL  =  SC  x  NPREDF  x  AMAX  x 


1-exp 


-EE  x  I  A  BSj  x  3600 
AMAX  x  SC 


(1.2) 

(1.3) 
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Table  2.  Parameters,  variables  and  constants,  grouped  according  to  model 
processes. 

Var/Constant 

c/va 

Value  SI 

Value  Sa 

Unit 

Description 

Source 

SI 

Source 

Sa 

Phenological  cycle  and  development 

FLV(T) 

v  (tab) 

0.190 

0.135 

Unitless 

Fraction  of  total  dry  matter  increase 
allocated  to  leaves  as  function  of  DVS 

1 

2,3,4 

FST(T) 

v  (tab) 

0.520 

0.421 

Unitless 

Fraction  of  total  dry  matter  increase 
allocated  to  stems  as  function  of  DVS 

1 

2,3,4 

FRT(T) 

v  (tab) 

0.290 

0.444 

Unitless 

Fraction  of  total  dry  matter  increase 
allocated  to  roots  as  function  of  DVS 

1 

2,  3,  3 

DDTMP 

V 

°C 

Daily  average  temperature  (field  site) 

DVRV(T) 

v  (tab) 

0.015 

0.022 

d'1 

Development  rate  after  flowering  as 
function  of  temperature 

Calibr. 

Calibr. 

DVRR(T) 

v  (tab) 

0.040 

0.015 

d'1 

DVR  prior  to  flowering  as  function  of 
temperature 

Calibr. 

Calibr. 

DVS 

V 

Unitless 

Development  phase 

Calibr. 

Cal;ibr. 

Wintering,  sprouting,  sprout  elongation 

NPL 

c 

30 

77 

m~2 

Plant  density 

5 

6 

NDTUB 

V 

101 

NA 

rrf2 

Dormant  tuber  density 

7 

INTUB 

c 

0.765 

NA 

g  DW  tuber 1 

Tuber  size 

8 

RDTU 

c 

0.014 

NA 

d'1 

Relative  tuber  death  rate  (on  number 
basis) 

7 

NTUBD 

V 

N  m'2 

Dead  tuber  number 

NTUBPD 

V 

N  m'2 

Dead  tuber  number  previous  day 

NGTUB 

V 

N  m'2 

Sprouting  tuber  number 

REMOB 

V 

g  CH20  rrf2  d'1 

Remobilization  rate  of  carbohydrates 

ROC 

0.0576 

0.0576 

g  CH20  g  1  DW  d'1 

Relative  conversion  rate  of 
tuber/rhizome  into  plant  material 

9 

9 

RCSHST 

c 

12 

NA 

m  g'1  DW 

Relation  coefficient  tuber  weight-stem 
length 

9 

CRIFAC 

c 

0.87 

0.535 

ci  DW  layer1  plant 

Critical  shoot  weight  per  0. 1  -m  depth 
layer 

1 

10 

SURPER 

c 

99 

NA 

d 

Survival  period  for  sprouts  without  net 
photosynthesis 

11 

TWGTUB 

V 

g  DW  rrf2 

Total  dry  weight  of  sprouting  tubers 

IWGRIZ 

c 

40.1 

928 

g  DW  rtf2 

Initial  rhizome  weight 

RDRIZ 

V 

0.00042 

0.00042 

d'1 

Relative  rhizome  death  rate 

CRRIZ 

c 

10.6 

400 

g  DW  rrf2 

Critical  rhizome  weight 

Photosynthesis,  maintenance,  growth,  and  assimilate  partitioning 

SC 

c 

J  rrf2  s'1 

Solar  constant  corrected  for  varying 
distance  sun-earth 

12 

12 

TL 

c 

0.1 

0.1 

m 

Thickness  depth  layer 

IABS(i) 

V 

J  m'2  s'1 

Total  irradiance  absorbed  by  depth 
layer  i 

IABSL(i) 

V 

J  m'2  s  1 

Total  irradiance  absorbed  by  shoots 
in  depth  layer  i 

IRZ(i) 

V 

J  rrf2  s'1 

Total  photosynthetically  active  part  of 
irradiance  on  top  of  depth  layer  i 

SCO) 

V 

g  DW  rrf2 

Shoot  dry  matter  in  depth  layer  i 

K(T) 

v  (tab) 

0.019 

0.00241 

m2  g  1  DW 

Plant  species  specific  light  extinction 
coefficient  as  function  of  DVS 

1 

10 

EPISHD 

V  (tab) 

0-0.43 

0-0.43 

Unitless 

Fraction  of  irradiation  shaded  by 
epiphytes 

13 

13 

AMX 

c 

0.0165 

0.0165 

g  C02  g  1  DW  h  1 

Potential  C02  assimilation  rate  at  light 
saturation  for  under-water  shoots 

14 

14 

(Continued) 
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Table  2.  Continued 

Var/Constant 

c/va 

Value  SI 

Value  Sa 

Unit 

Description 

Source 

SI 

Source 

Sa 

AMX2 

c 

0.0620 

0.011 

g  C02  g  1  DW  h’1 

Potential  C02  assimilation  rate  at  light 
saturation  for  above-water  shoots 

1,  15 

4,  15 

AMAX 

V 

g  C02  g  1  DW  h'1 

Actual  C02  assimilation  rate  at  light 
saturation  for  shoots 

EE 

c 

0.000011 

0.000014 

g  co2  J  1 

Initial  light  use  efficiency  for  shoots 

12 

12 

NPRAT 

v  (tab) 

6-8 

5-8 

Unitless 

Plant  biomass  N:P  ratio 

NPREDF(T) 

v(eq) 

0-1 

0-1 

Unitless 

Relative  AMX  factor  to  account  for 
nutrient  limitation 

REDF(T) 

v  (tab) 

1 

1 

Unitless 

Relative  reduction  factor  for  AMX  to 
account  for  senescence  plant  parts 

User 

def. 

User 

def. 

REDAM 

c 

1 

1 

Unitless 

Relative  reduction  factor  to  relate 

AMX  to  water  pH  and  oxygen  level 

User 

def. 

User 

def. 

REDAM1 

v  (tab) 

0-1 

0-1 

Unitless 

Relative  reduction  factor  to  relate 

AMX  to  water  current  velocity 

AMTMP(T) 

v  (tab) 

0-1 

0-1 

Unitless 

Daytime  temperature  effect  on  AMX 
as  function  of  DVS 

16 

16 

FGL 

V 

g  C02  rrf2  h'1 

Instantaneous  C02  assimilation  rate 
per  vegetation  layer 

GPHOT 

V 

g  CH20  m’2  d’1 

Daily  total  gross  assimilation  rate  of 
the  vegetation 

DMPC(T) 

v  (tab) 

0-1 

0-1 

Unitless 

Dry  matter  allocation  to  each  plant 
layer 

1 

10 

ASRQ 

V 

g  CH20  g  1  DW  d’1 

Assimilation  requirement  for  plant  dry 
matter  production 

FL(T) 

v  (tab) 

0-1 

0-1 

Unitless 

Leaf  dry  matter  allocation  to  each 
layer  of  shoot  as  function  of  DVS 

1 

2,3,4 

GLV 

V 

g  DW  m’2  d’1 

Dry  matter  growth  rate  of  leaves 

GST 

V 

g  DW  m’2  d’1 

Dry  matter  growth  rate  of  stems 

GRT 

V 

g  DW  m’2  d’1 

Dry  matter  growth  rate  of  roots 

GTW 

V 

g  DW  m’2  d’1 

Dry  matter  growth  rate  of  the  vegeta¬ 
tion  (excl.  tubers,  rhizomes) 

TWLVG 

V 

g  DW  m’2 

Total  dry  weight  live  leaves 

TWSTG 

V 

g  DW  m’2 

Total  dry  weight  live  stems 

TWRTG 

V 

g  DW  m’2 

Total  dry  weight  live  roots 

TGW 

V 

g  DW  m’2 

Total  live  plant  dry  weight  (excl. 
tubers,  rhizomes) 

MAINT 

V 

g  CH20  m’2  d’1 

Maintenance  respiration  rate  vegeta¬ 
tion 

MAINTS 

V 

g  CH20  m’2  d'1 

Maintenance  respiration  rate  vegeta¬ 
tion  at  reference  temperature 

Upper  biomass 
limit 

c 

950 

2,250 

g  DW  m*2 

Maximum  plant  biomass 

7 

2,  17 

Flowering,  translocation,  senescence,  and  formation  of  wintering  organs 

RTR 

c 

0.150 

NA 

g  DW  tuber 1  d'1 

Maximum  relative  tuber  qrowth  rate  at 
20°C 

1 

RTRL 

V 

NA 

g  DW  tuber’1  d'1 

Relative  tuber  growth  rate  at  ambient 
temperature 

CVT 

c 

1.05 

1.05 

Unitless 

Conversion  factor  for  translocated  dry 
matter  into  CFI20 

12 

12 

NINTUB 

c 

2.0 

NA 

N  plant’1 

Tuber  number  concurrently  initiated 
per  plant 

1 

TWCTUB 

c 

45.9 

NA 

g  DW  m’2 

Total  critical  dry  weight  of  new  tubers 

1,  5,  8 

NNTUB 

V 

NA 

N  m’2 

New  tuber  number 

RDR(T) 

v  (tab) 

0.030 

0.030 

d’1 

Relative  death  rate  of  leaves  as  func¬ 
tion  of  DAVTMP  (on  DW  basis) 

1 

6 

(Continued) 
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Table  2.  Concluded 

Var/Constant 

c/va 

Value  SI 

Value  Sa 

Unit 

Description 

Source 

SI 

Source 

Sa 

RDS(T) 

v  (tab) 

0.030 

0.030 

Q. 

Relative  death  rate  of  stems  and  roots 
as  function  of  DAVTMP  (on  DW 
basis) 

1 

6 

TEFF(T) 

v  (tab) 

Unitless 

Relative  effective  temperature 
function  accounting  for  temperature 
effect  on  maintenance  respiration, 
remobilization,  maximum  tuber  growth 
and  death  rates  as  function  of 
temperature 

Calibr. 

TRANS 

V 

g  CH20  m'2  d'1 

Translocation  rate  of  carbohydrates 

I .  K.P.  Kenow,  unpublished  results  2008;  2.  Gallagher  et  al.  1984;  3.  Smart  and  Barko  1980;  4.  Best  et  al.  2008b;  5.  Low  and 
Bellrose  1944;  6.  Dai  and  Wiegert  1996b;  7.  Clark  and  Clay  1985;  8.  Marburger  1993;  9.  Bowes  et  al.  1979;  10.  Morris  1989; 

II.  Best  and  Boyd  2001;  12.  Penning  de  Vries  and  Van  Laar  1982a,b;  13.  Best  et  al.  2005;  14.  Titus  and  Adams  1979; 

15.  Penning  de  Vries  and  Van  Laar  1982a,  b;  16.  Best  and  Boyd  2003;  17.  Dai  and  Wiegert  1996a. 

“Ac  indicates  that  the  parameter  is  a  constant.  A  v  indicates  a  variable,  eq  and  tab  indicate  that  the  parameter  is  implemented  in 
the  model  as  an  equation  and  a  table,  respectively.  Abbreviations:  SI  =  S.  latifolia ;  Sa  =  S.  alterniflora. 

METHODS  SIMULATION  STUDIES:  The  generic  model  is  composed  by  a  framework,  in 
which  the  model  applications  ARROW  and  CORDG  can  be  run  simultaneously.  The  model 
requires  daily  values  of  the  following  environmental  variables  as  inputs:  water  depth,  water 
transparency,  temperature  (water  or  air)  and  irradiance.  Among  the  required  inputs,  the  data  on 
water  depth  can  be  derived  from  local  and  regional  stage  observations  obtained  from  a  web- 
based  database.  Data  on  water  transparency  can  be  derived  from  Secchi  disk  observations  also 
obtained  from  a  web-based  database  using  the  relationship  of  Giessen  et  al  (1990).  Thus,  the 
light  extinction  coefficient  (L),  required  as  input  for  these  ecological  models,  can  be  derived 
from  measured  Secchi  disk  depths  following  L  (in'1)  =  1.65/  Secchi  disk  depth  (m).  The  latter 
relationship  is  valid  for  turbid,  shallow  water  only.  Both  water  depth  and  water  transparency  data 
can  also  be  derived  from  hydrodynamic  and  sediment  transport  model  results  (Best  et  al.  2008a). 
Data  on  irradiance  and  air  temperature  can  be  obtained  from  local  or  regional  weather  stations. 

In  the  present  simulation  studies  the  following  environmental  data  were  used  as  inputs:  (i)  a  con¬ 
stant  water  depth  of  0.2  m  unless  indicated  otherwise  (Table  3;  0.2  m  depth  is  typical  for  shallow 
water  bodies  such  as  river  pools,  but  usually  daily  and  seasonal  fluctuations  occur,  as  docu¬ 
mented  by  Best  and  Boyd  2008);  (ii)  light  extinction  coefficients,  typical  for  turbid  water  such  as 
river  pools  and  peat  lakes  (1.81-2.0  m  ,  Table  3);  (iii)  weather  data,  either  typical  for  a  temper¬ 
ate  climate  (Mississippi  River  Pool  8  and  Pool  9,  near  La  Crosse,  WI,  for  years  1982  (calibra¬ 
tion)  and  2006  (validation)  in  which  field  data  on  S.  latifolia  biomass  distribution  were  collected; 
Sapelo  Island,  GA,  for  years  1991  (calibration)  and  1972  (validation)  in  which  field  data  on 
S.  alterniflora  biomass  distribution  were  collected),  or  typical  for  a  near  subtropical  climate 
where  both  species  also  abundantly  grow  but  descriptions  of  wax,  wane,  and  coexistence  are  still 
lacking,  i.e.,  Kenner  (near  New  Orleans),  Louisiana,  2006  (exploration;  Table  3). 
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Table  3.  Variables  and  constants,  grouped  according  to  field  site  characteristics 

and  management. 

Var/Constant 

c/va 

Value  SI 

Value 

Sa 

Unit 

Description 

Source  SI 

Source  Sa 

Field  site  characteristics'1 

DPT(T) 

v  (tab) 

0.07-2.09 

0.2 

m 

Water  depth  (field  site) 

1,  User 
def. 

User  def. 

WTMP(T) 

v  (tab) 

°c 

Daily  water  temperature  as  function  of 
day  no  (field  site) 

User  def. 

User  def. 

L(T) 

V(tab) 

1.81 

2.0 

m1 

Water  type  specific  light  extinction 
coefficient  as  function  of  day  no  (field 
site) 

1,  User 
def. 

User  def. 

WVEL 

v  (tab) 

0-100 

0-100 

cm  s'1 

Water  type  specific  current  velocity  as 
function  of  day  no  (field  site) 

User  def. 

User  def. 

TGWM(T) 

v  (tab) 

g  DW  m 2 

Total  live  dry  weight  measured  as 
function  of  day  no  (field  site) 

User  def. 

2,  3,  User 
def. 

NTM(T)  -Va 

v  (tab) 

101 

NA 

N  rrT2 

Tuber  density  measured  as  function  of 
day  no  (field  site) 

User  def. 

Management  (harvesting) 

HAR 

c 

0  or  1 

0  or  1 

Harvesting  switch  (0=off,  1=on) 

User  def. 

User  def. 

HARDAY 

c 

1-365 

1-365 

d 

Harvesting  day  number 

User  def. 

User  def. 

HARDEP 

c 

m 

Harvesting  depth  (measured  in  0.1 -m 
increments  from  water  surface) 

User  def. 

User  def. 

1 .  Clark  and  Clay  1 985;  2.  Gallagher  et  al.  1 984;  3.  Dai  and  Wiegert  1 996a. 

a  A  c  indicates  that  the  parameter  is  a  constant.  A  v  indicates  that  a  variable  and  tab  indicates  that  the  parameter  is  implemented 
in  the  model  as  a  table.  Abbreviations:  SI  =S.  latifolia-,  Sa=  S.  alterniflora. 

b  Temperate  field  site:  For  SI,  La  Crosse,  Wisconsin  (lat  43°  10’N,  long  91°  30’W);  weather  file  1982;  For  Sa,  Brunswick,  Georgia 
(lat  31  °  1 5’N,  long  81  °  28’W);  weather  file  1 991 ;  Near  subtropical  field  site:  Kenner,  Louisiana  (lat  29°  59’N,  long  90°  1 5’W); 
weather  file  2006. 

SIMULATIONS 

S.  latifolia  Base  Runs.  ARROW  was  selected  from  the  framework  and  run  using  the  nominal 
parameter  values  (Table  2),  field  site  variables  (Table  3)  and  weather  data  (irradiance  and  air 
temperature,  La  Crosse,  WI;  lat  43°  30’N,  long  91°  10’W)  as  inputs  for  a  1-year  period.  The 
simulated  biomass  of  S.  latifolia  plants  (composed  by  shoots  and  roots),  rhizomes  and  tubers  are 
shown  in  Figure  3.  Simulated  plant  biomass  augmented  with  rhizome  biomass  compared  well 
with  the  measured  ‘plant’  biomass,  in  which  plants  and  rhizomes  had  not  been  separated  during 
harvesting.  Simulated  plant  biomass  reached  its  maximum  100  days  earlier  than  rhizome  bio¬ 
mass,  and  measured  “plant”  biomass  reached  its  maximum  10  days  after  simulated  plant  bio¬ 
mass.  A  comparison  of  simulated  and  measured  rhizome  biomass  was  not  possible  because  rhi¬ 
zome  biomass  had  not  been  determined  in  the  calibration  data  set.  Simulated  tuber  number 
decreased  during  summer  and  increased  from  August  onwards  to  a  level  that  was  somewhat 
greater  than  initially.  The  thus-calibrated  model  simulated  a  stable  S.  latifolia  population  that 
persisted  by  sprouting  from  rhizomes  as  well  as  from  tubers  in  spring,  but  did  not  attain 
280  tubers  m'2  as  reported  by  Clark  and  Clay  (1985)  in  Pool  9  in  1982. 


12 


ERDC  TN-SWWRP-1 1-1 
January  201 1 


Figure  3.  Simulated  biomass  of  plants,  rhizomes,  and  tubers  (A),  tuber  numbers  (B),  measured  plants 
(C),  and  measured  tuber  numbers  (D)  of  Sagittaria  latifolia  in  Upper  Mississippi  Pool  9,  IA. 
Nominal  run.  Field  data  1982  from  Clark  and  Clay  (1985);  climatological  data  1982, 

La  Crosse,  Wisconsin  (lat  43°  30’N,  long  91°  10’W);  water  depth  0.1 1  to  2.10  m;  light 
extinction  coefficient  1.81. 


S.  latifolia  Validation  Runs.  Results  of  the  validation  run  indicated  that  simulated  plant  bio¬ 
mass  was  similar  to  measured  plant  biomass,  but  that  in  this  case  simulated  tuber  numbers 
greatly  exceeded  the  measured  ones  (Figure  4). 


Figure  4.  Simulated  biomass  of  plants,  rhizomes,  and  tubers  (A),  tuber  numbers  (B),  measured  plants 
(C),  and  measured  tuber  numbers  (D)  of  Sagittaria  latifolia  in  Upper  Mississippi  Pool  8,  Wl. 
Validation  run.  Field  data  2006  from  K.  P.  Kenow  (unpublished  Long  Term  Research 
Management  Program,  2006);  climatological  data  2006,  La  Crosse,  Wisconsin  (lat  43°  30’N, 
long  91  °  1 0’W);  water  depth  0  to  0.62  m;  light  extinction  coefficient  0.98  to  2.54. 
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The  discrepancies  between  the  simulated  and  the  measured  tuber  numbers  in  the  calibration  and 
validation  runs  could  be  explained  in  two  different  ways.  On  the  one  hand,  the  tuber  data  were 
extremely  scarce  and  representative  for  only  one  or  two  points  in  time  during  the  year  (i.e.,  on 
30  June  and  15  September  in  1982,  and  on  30  June  in  2006),  making  a  positive  correlation 
unlikely.  On  the  other  hand,  since  measured  tuber  density  was  high  at  the  end  of  summer  in  Pool 
9  in  1982  and  low  in  Pool  8  in  2006  (the  tubers  were  measured  at  sites  within  waterfowl 
exclosures  in  which  grazing  was  prevented),  it  could  not  be  ruled  out  that  plant  populations 
produce  tuber  numbers  which  differ  greatly  between  years  in  a  temperate  climate,  making 
rhizomes  extremely  important  organs  for  population  persistence  at  this  latitude.  The  latter 
possibility  was  further  explored  by  conducting  model  runs  for  four  different  water  level  -  year 
combinations  (Figure  5).  The  results  of  these  simulations  indicate  that  simulated  plant-rhizome 
biomass  varied  by  a  factor  of  1.5,  and  tuber  numbers  were  sufficient  in  all  cases  to  enable 
persistence  of  the  population  and  sprouting  from  both  rhizomes  and  tubers,  rendering  a  large 
variation  in  tuber  numbers  between  years  unlikely. 


Figure  5.  Simulated  biomass  of  plants  (including  rhizomes)  and  tubers  of  Sagittaria  latifolia  (upper)  and 
measured  typical  water  level  fluctuations  (lower)  in  Upper  Mississippi  System  Pool  8, 
Wisconsin.  Nominal  initial  biomass  and  light  extinction  coefficient  values;  climatological  data 
La  Crosse,  Wisconsin  (lat  43°  30’N,  long  91°  10’W). 
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S.  latifolia  Runs  in  a  Subtropical  Climate.  To  investigate  whether  ARROW  could  be  used 
to  simulate  behavior  of  a  S.  latifolia  community  in  a  subtropical  climate,  a  run  was  conducted  for 
a  more  southern  site,  Kenner,  Louisiana  (lat  29°  59’N,  long  90°  15’W)  (Figure  6).  Results  of  this 
run  indicated  that  simulated  plant  biomass  exceeded  plant  biomass  in  a  temperate  climate  by  a 
factor  of  four,  sufficient  rhizome  biomass  and  large  numbers  of  tubers  were  produced.  From  this 
information,  it  can  also  be  concluded  that  the  model  simulated  a  stable  S.  latifolia  population  that 
persisted  by  sprouting  from  rhizomes  and  from  tubers  in  spring. 


Figure  6.  Simulated  biomass  of  plants,  rhizomes,  and  tubers  (A),  and  tuber  numbers  (B),  of  Sagittaria 
latifolia  in  a  more  southern  climate.  Nominal  initial  biomass,  light  extinction  coefficient,  and 
water  depth  values;  climatological  data  Kenner,  Louisiana  (lat  29°  59’N,  long  90°  15’W).  For 
comparison,  measured  biomass  of  plants  (C),  and  tuber  numbers  (D)  in  Upper  Mississippi 
Pool  9,  IA  (near  La  Crosse,  Wisconsin;  lat  43°  30’N,  long  91°  10’W)  provided. 

S.  alterniflora  Base  Runs.  CORDG  was  selected  from  the  framework  and  run  using  the 
nominal  parameter  values  pertaining  to  a  S.  alterniflora  population  in  Sapelo  Island,  Georgia 
(Table  2).  Field  site  variables  (Table  3)  and  weather  data  (irradiance  and  air  temperature,  Bruns¬ 
wick,  Georgia;  lat  31°  15’N,  long  81°  28’ W)  were  used  as  inputs  for  a  1-year  period.  CORDG 
was  calibrated  after  the  ‘tall’  S.  alterniflora  variety,  which  usually  grows  relatively  close  to  mean 
high  water  level  and  along  creeks.  The  ‘tali’  variety  differs  in  that  it  has  a  greater  plant  species 
characteristic  light  extinction  coefficient  for  shoots  (0.00241  m  g'  DW;  Morris  1989)  than  the 
‘short’  variety  (0.00187  nr  g'1  DW;  Morris  1989;  Dai  and  Wiegert  1996a),  which  grows  at  a 
higher  elevation  of  the  marsh.  Simulated  biomass  of  S.  alterniflora  plants  (including  shoots  and 
roots),  roots,  and  rhizomes  is  shown  in  Figure  7.  Simulated  biomass  showed  more  variation  with 
season  than  measured  biomass,  and  there  was  an  overall  good  agreement  between  both.  Simu¬ 
lated  plant  biomass  reached  two  maxima,  one  at  the  end  of  March  and  one  at  the  end  of  October. 
Measured  plant  (shoot  plus  root)  biomass  showed  one  maximum  at  the  end  of  August,  was  only 
determined  at  three  points  in  time  (in  January,  end  of  August  and  end  of  December),  and  mean 
values  ±  standard  deviations  matched  simulated  plant  biomass  values.  Measured  data  on  roots 
and  rhizomes  were  scarce,  without  replicates,  inhibiting  unequivocal  matching  of  measured  and 
simulated  values  -  but  measured  and  simulated  biomass  were  in  the  same  order  of  magnitude  as 
were  the  trends.  The  thus-calibrated  model  simulated  a  stable  S.  alterniflora  population  that 
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persisted  by  sprouting  from  rhizomes  in  spring,  with  aboveground  biomass  (composed  by 
shoots)  reaching  a  lower  maximum  of  900  g  DW  irf2  than  1500  g  DW  m"2  of  belowground 
biomass  (composed  by  rhizomes  and  roots),  and  rhizome  mass  exceeding  plant  mass  in  winter. 


S.  alterniflora 


Time  (d) 


Time  (d) 


Time  (d) 


Figure  7.  Simulated  biomass  of  plants,  roots  and  rhizomes  (A),  rhizomes  (B),  and  measured  shoots, 
roots  and  rhizomes  (C)  of  Spartina  alterniflora  on  Sapelo  Island,  Georgia.  Nominal  run.  Field 
data  1991  from  Dai  and  Wiegert  (1996  a,  b);  climatological  data  1991,  Brunswick,  Georgia 
(lat  31°  15’N,  long  81°  28’W);  water  depth  0.20  m;  light  extinction  coefficient  2.0. 


S.  alterniflora  Validation  Runs.  Results  of  the  validation  run  indicated  that  simulated  shoot 
biomass  (plants  minus  roots)  was  similar  to  measured  shoot  biomass  (Figure  8).  It  was  not  possi¬ 
ble  to  compare  simulated  to  measured  belowground  biomass,  because  no  belowground  biomass 
was  determined  in  the  validation  data  set. 
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Figure  8.  Simulated  biomass  of  plants,  roots  and  rhizomes  (A),  rhizomes  (B),  and  measured  shoots 
(C)  of  Spartina  altemiflora  on  Sapelo  Island,  Georgia.  Validation  run.  Field  data  1972  from 
Gallagher  etal.  (1980);  climatological  data  1972,  Brunswick,  Georgia  (Iat31°  15’N,  long  81° 
28’W);  water  depth  0.20  m;  light  extinction  coefficient  2.0. 


S.  altemiflora  Runs  in  a  Subtropical  Climate.  To  investigate  whether  CORDG  could  be 
used  to  simulate  behavior  of  a  S.  altemiflora  community  in  a  subtropical  climate,  a  run  was  con¬ 
ducted  for  the  same,  more  southern  site  as  was  used  to  test  ARROW,  (i.e.,  Kenner,  Louisiana) 
(Figure  9).  Results  of  this  run  indicated  that  simulated  plant,  root,  and  rhizome  biomass  in 
Kenner  were  a  factor  of  1.4  greater  than  those  on  Sapelo  Island:  the  model  simulated  a  stable  S. 
altemiflora  population  that  persisted  by  sprouting  from  rhizomes  in  spring.  Results  of  an 
additional  run  for  the  “short”  plant  variety,  which  forms  a  large  part  of  CL  marshes,  indicated  an 
1 8  percent  reduction  in  biomass  production  compared  to  production  by  the  “tall”  variety  by  a 
stable  plant  population. 

Sensitivity  Analysis.  Sensitivity  analysis  of  a  simulation  model  must  be  undertaken  to  assess 
the  parameters  likely  to  strongly  affect  model  behavior.  The  current  analysis  is  based  on  the 
effect  of  a  change  in  one  parameter  while  all  other  parameters  are  kept  the  same.  The  parameter 
under  study  was  changed  and  1-year  simulations  were  conducted  under  nominal  environmental 
conditions.  The  nominal  parameter  values  (as  presented  in  Table  2)  were  chosen  as  a  reference 
level.  The  results  were  compared  with  those  of  a  nominal  run.  Each  parameter  was  increased 
once  by  20  percent  and  decreased  once  by  20  percent.  The  relative  sensitivity  (RS)  of  a  parame¬ 
ter  was  then  defined  as  the  relative  change  in  the  variable  on  which  the  effect  was  tested  divided 
by  the  relative  change  in  the  parameter  (Ng  and  Loomis  1984).  The  effects  of  thirteen  parameters 
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Figure  9.  Simulated  biomass  of  plants,  roots  and  rhizomes  (A  =  tall,  C  =  short  vegetation),  and  rhizomes 
(B  =  tall,  D  =  short  vegetation),  of  Spartina  alterniflora  in  a  more  southern  climate.  Nominal 
initial  biomass,  light  extinction  coefficient,  and  water  depth  values;  climatological  data  Kenner, 
Louisiana  (lat  29°  59’N,  long  90°  15’W). 


on  three  state  variables,  representing  different  plant  biomass  compartments,  were  tested.  A 
model  variable  is  considered  sensitive  to  a  change  in  the  value  of  a  parameter  at  RS  >  0.5  and 


<-0.5. 


( yield i  -  yield,.) 


RS  = 


yield,. 


( parami  -  paramr ) 

param,. 


where 

yield,  =  value  at  parameter  value  i 
yield,-  =  value  at  reference  parameter  value  r 
parami  and  param,-  =  as  above 

The  three  state  variables,  maximum  plant  biomass,  end-of-year  tuber  number  and  end-of-year 
rhizome  biomass  in  the  model  were  sensitive  to  parameter  changes,  particularly  those  affecting 
carbon  capture  (AMX2- potential  CO2  assimilation  rate  at  light  saturation  for  above-water  shoots; 
and  EE-  initial  light  use  efficiency  for  shoots;  Table  4).  In  addition,  end-of-year  rhizome  and/or 
tuber  number  was  sensitive  to  development  rates  ( DVRVT ,  DVRRT)  and  relative  death  rate 
( RDR )  and  in  SLAT  only  to  plant  density  ( NPL ),  initial  rhizome  weight  ( IWGRIZ)  and  relative 
tuber  growth  rate  ( RTR ). 
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Table  4.  Relative  sensitivity  (RS)  of  selected  state  variables  to  deviations  in 
parameter  values  from  their  nominal  values  as  presented  in  Table  2.  The  RS  of  a 
parameter  is  the  relative  change  in  the  variable  on  which  the  effect  was  tested 
divided  by  the  relative  change  in  the  parameter.  A  model  variable  is  considered 
sensitive  to  a  change  in  the  value  of  a  parameter  at  RS  >  0.5  and  <  -0.5.  Results 
were  obtained  in  1 -year  simulations  under  nominal  conditions.  Sensitive  values 
shaded. 

Parameter 

Relative  sensitivity 

Name 

Value 

S.  latifolia 

S.  alterniflora  j 

Max.  plant 
biomass 

EOY1  rhiz. 
biomass 

EOY1  tuber 
no 

Max.  plant 
biomass 

EOY1  rhiz. 
biomass 

DVRVT 

+20% 

0.25 

0.48 

0.12 

-0.11 

-0.09 

-20% 

0.42 

4.99 

2.56 

-0.07 

-0.15 

DVRRT 

+20% 

0.25 

0.48 

0.12 

-0.68 

-0.76 

-20% 

0.24 

2.85 

1.46 

-1.26 

-0.19 

NPL 

+20% 

0.04 

0.01 

0.85 

-0.55 

-0.38 

-20% 

0.04 

0.01 

-0.90 

-0.04 

-0.03 

INTUB 

+20% 

0.04 

0.01 

0.0 

NA 

NA 

-20% 

0.04 

0.01 

-2.20 

NA 

NA 

ROC 

+20% 

0.03 

0.00 

-0.01 

0.05 

0.04 

-20% 

0.03 

0.01 

-0.02 

0.00 

0.00 

RTR' 

+20% 

-0.07 

-4.99 

-0.98 

NA 

NA 

-20% 

-0.07 

-0.07 

0.0 

NA 

NA 

RDTU 

+20% 

0.00 

0.00 

-0.14 

NA 

NA 

-20% 

0.00 

0.00 

-0.20 

NA 

NA 

IWGRIZ 

+20% 

0.12 

0.02 

0.00 

0.20 

0.20 

-20% 

0.17 

4.99 

2.49 

0.02 

0.10 

TRAFAC 

+20% 

-0.30 

0.14 

0.02 

-0.13 

0.37 

-20% 

-0.30 

4.99 

2.54 

-0.33 

0.31 

A  MX 

+20% 

0.00 

0.00 

0.00 

0.00 

0.00 

-20% 

0.00 

0.00 

0.00 

0.00 

0.00 

AMX2 

+20% 

1.11 

0.00 

0.00 

0.75 

0.78 

-20% 

0.69 

4.99 

2.34 

2.05 

2.69 

EE 

+20% 

0.90 

0.92 

2.49 

1.46 

1.42  i 

-20% 

1.02 

3.66 

1.35 

2.04 

2.75 

RDR 

+20% 

-0.17 

-4.99 

-2.48 

-1.95 

-2.28 

-20% 

-0.21 

-0.53 

0.00 

-0.41 

-0.98 

Note:  1  EOY  -  end  of  year 

The  sensitivity  of  maximum  plant  biomass  and  end-of-year  tuber  number/rhizome  biomass  to 
changes  in  environmental  factors  was  assessed  by  following  the  same  approach  as  for  sensitivity 
analysis  of  the  model  parameters.  For  this  purpose,  parameter  changes  were  based  on  value 
ranges  taken  from  literature,  which  sometimes  differed  more  than  20  percent  from  the  nominal 
parameter  values  presented  in  Table  4.  Also  in  this  analysis,  the  model  proved  to  be  sensitive 
(Table  5).  All  state  variables  were  sensitive  to  changes  in  climate.  Maximum  plant  biomass  was 
less  sensitive  in  SLAT  than  in  SPALT.  Changes  in  the  light  reflection  coefficient  at  the  water 
surface  ( RC)  had  no  significant  effect,  while  changes  in  the  other  tested  parameter  values  had 
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effects  which  decreased  in  the  order  of  the  water  type’s  specific  light  extinction  coefficient  (LT) 
greater  than  water  depth  ( DPTT ). 


Table  5.  Environmental  factor  analysis,  expressed  as  relative  sensitivity  (RS)  of 
selected  state  variables  to  deviations  in  parameter  values  from  their  nominal 
values,  as  presented  in  Table  2.  The  RS  of  a  parameter  is  the  relative  change  in 
the  variable  on  which  the  effect  was  tested  divided  by  the  relative  change  in  the 
parameter.  A  model  variable  is  considered  sensitive  to  a  change  in  the  value  of  a 
parameter  at  RS  >  0.5  and  <  -0.5.  Results  were  obtained  in  1-year  simulations 
under  nominal  conditions.  Sensitive  values  shaded. 

Parameter 

Relative  sensitivity 

S.  latifolia 

S.  alterniflora  j 

Max.  plant 

EOY3  rhiz. 

EOY3  tuber 

Max.  plant 

EOY5  rhiz.  j 

Name 

Value 

biomass 

biomass 

no 

biomass 

biomass 

Climate1 

Lat  43°30’  N 

Lat  29°59’  N 

-11.33 

-3.31 

1.35 

Lat  31  °N 

Lat  29°59’  N 

-2.05 

-1.58 

RC4 

1.00 

-0.06 

-0.06 

-0.06 

-0.03 

-0.03 

0.00 

-0.06 

-0.05 

0.00 

-0.11 

-0.10 

LT 

+20% 

-0.44 

-4.99 

-2.47 

-1.13 

-0.93 

-20% 

-0.45 

-0.28 

0.00 

-0.62 

-0.62 

DPTT 

+20% 

-0.15 

-4.99 

-2.49 

0.00 

0.00 

-20% 

-0.23 

-0.03 

0.00 

0.00 

0.00 

1 1  Climates  at  La  Crosse,  Wisconsin  (lat  43°  30'  N),  Kenner,  Louisiana  (lat  29°  59’  N) 

Brunswick,  Georgia  (lat.  31°  N) 

EOY  -  end  of  year 

4  Light  reflection  coefficient  at  water  surface;  to  enable  calculation  of  the  relative  sensitivity,  a  very  low  value  of  0.000001  was 

used. 

SUMMARY:  A  dynamic  simulation  modelling  approach  to  emergent  plant  biomass  fonnation 
has  been  developed  to  provide  a  tool  for  water  resource  managers.  It  is  now  possible  to  evaluate 
key  environmental  conditions  in  which  emergent  aquatic  vegetation  would  persist  or  produce 
excessive  biomass  with  ensuing  consequences  for  the  systems  in  which  they  grow,  whether  they 
are  affected  by  management  measures  or  not.  The  generic  model  is  composed  a  framework,  in 
which  the  model  applications  ARROW  for  Sagittaria  latifolia  and  CORDG  for  Spartina  alter- 
niflora  can  be  run  simultaneously.  The  model  describes  major,  carbon  flow-based  ecophysio- 
logical  processes  and  biomass  dynamics  of  two  common  plant  species.  It  also  contains  unique 
descriptions  of:  (1)  species-characteristic  vertical  distribution  of  shoot  biomass  which  enables  the 
calculation  of  the  fraction  of  irradiance  actually  available  for  absorption  by  the  plant;  (2)  recal¬ 
culation  procedures  of  this  vertical  distribution  with  daily  changes  in  water  level  and/or  shoot 
mass  removal  at  various  heights  and  levels  within  the  water  column,  which  enables  the  evalua¬ 
tion  of  regrowth  potential;  and  (3)  relationships  of  plant  process  parameters  with  site-specific 
climate  which  enables  the  evaluation  of  effects  of  different  climates.  Generally,  a  good  fit  was 
found  between  simulated  and  measured  biomass  in  the  field.  Sensitivity  analysis  showed  that  the 
model  is  very  sensitive  to  changes  in  process  parameters  influencing  carbon  flow. 
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PRODUCT  DEVELOPMENT  AND  AVAILABILITY:  Aquatic  plant  growth  models  are  avail¬ 
able  to  both  U.S.  Army  Corps  of  Engineers  (USACE)  and  non-USACE  interested  parties.  The 
model  can  be  downloaded  from  the  following  URL:  http: //el. erdc. usace. army. mil/products. 
cfm?Topic=model&Type=aauatic.  Model  descriptions  and  user  manuals  can  be  downloaded 
from  the  same  web-page. 

POINTS  OF  CONTACT:  This  technical  note  was  prepared  by  Dr.  Elly  P.H.  Best  and 
William  A.  Boyd,  research  biologist  and  mathematician,  respectively,  U.S.  Army  Engineer 
Research  and  Development  Center,  Environmental  Laboratory;  and  Kevin  P.  Kenow,  Upper 
Midwest  Environmental  Sciences  Center,  USGS,  La  Crosse,  Wisconsin.  The  model  was 
conducted  as  an  activity  of  the  ecological  model  development  work  unit  of  the  System-Wide 
Water  Resources  Program  (SWWRP).  For  information  on  SWWRP,  please  visit 
h ftps ://swwrp. usace.armv.mil/  or  contact  the  Program  Manager,  Dr.  Steven  L.  Ashby  at 
Steven. L. Ashby  @  erdc.  us  ace,  army. mil.  Questions  about  this  technical  note  may  be  addressed  to 
Mr.  Boyd  at  (601-634-3705;  William. A.Boyd(/d), usace. army,  mil).  This  technical  note  should  be 
cited  as  follows: 

Best,  E.  P.  H.,  W.  A.  Boyd,  and  K.  P.  Kenow.  2011.  A  Generic  Modeling 
Approach  to  Biomass  Dynamics  of  Sagittaria  latifolia  and  Spartina  altemiflora. 
SWWRP  Technical  Notes  Collection.  ERDC  TN-SWWRP-1 1-1,  Vicksburg,  MS: 

U.S.  Army  Engineer  Research  and  Development  Center,  https ://swwrp. usace. 
army.mil/ 
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